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ABSTRACT 



Context. We present a self-consistent model of solar coronal heating in which we include the dynamical effect of the background 
magnetic field along a coronal structure by using exact results from wave MHD turbulence. 

Aims. We evaluate the heating rate and the microturbulent velocity for comparison with observations in the quiet corona, active re- 
gions and also coronal holes. 

Methods. The coronal structures are assumed to be in a turbulent state maintained by the slow erratic motion of the magnetic foot- 
points. A description of the large-scale and the unresolved small-scale dynamics are given separately. From the latter, we compute 
exactly (or numerically for coronal holes) turbulent viscosites used in the former to self-consistently close the system and derive the 
heating flux expression. 

Results. We show that the heating rate and the turbulent velocity compare favorably with coronal observations. 

Conclusions. Although the Alfven wave turbulence regime is strongly anisotropic, and could reduce a priori the heating efficiency, it 

provides a unexpected satisfactory model of coronal heating for both magnetic loops and open magnetic field lines. 
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1. Introduction 

Information about the solar corona from spacecraft missions like 
Yohkoh, SoHO (Solar & Heliospheric Observatory) or TRACE 
(Transition Region And Coronal Explorer) launched in the 
1990s reveals a very dynamical and complex medium structured 
in a network of magnetic field lines. These observations clearly 
demonstrated the fundamental role of the magnetic field on the 
plasma dynamics in the solar atmosphere. The solar corona con- 
tains a variety of structures over a broad range of scales from 
about 10 5 km until the limit of resolution (about one arcsec). It 
is very likely that structures at much smaller scales exist but have 
not yet been detected. The new spacecrafts STEREO, Hinode or 
SDO (Solar Dynamics Observatory) will help our understanding 
of the small-scale nature of the corona. 

Observations in UV and X-ray show a solar corona that is 
extremely hot with temperatures exceeding 10 6 K - close to hun- 
dred times the solar surface temperature. These coronal tempera- 
tures are highly inhomogeneous: in the quiet corona much of the 
plasma lies near 1-2 x 10 6 K and 1-8 x 10 6 K in active regions. 
Then, one of the major questions in solar physics concerns the 
origin of such high values of coronal temperature. The energy 
available in the photosphere is clearly sufficient to supply the to- 
tal coronal losses (Withbroe & Noyes 1977) which is estimated 
to be 10 4 J m _2 s _1 for active regions and about one or two orders 
of magnitude smaller for the quiet corona and coronal holes. The 
main issue is thus to understand how the available photospheric 
energy is transferred and accumulated in the solar corona, and 
by what processes it is dissipated. 

It is widely believed that the energy input comes from 
the slow random motion of the convective layer below the 
photosphere, but the mechanisms that heat the solar corona 



remain controversial. Heating models are often classified 
into two categories: "AC" and "DC" heating. "AC" heating 
by Alfven waves was suggested for the first time by Alfven 
(1947). In non-uniform plasmas, this heating mechanism could 
be sustained by the resonant absorption of Alfven waves 
(Hollweg 1984), or by phase mixing (Heyvaerts & Priest 
1983). The necessary condition for "AC" heating is that the 
characteristic time of motion excitation has to be shorter than 
the characteristic time of wave propagation across the coronal 
loops. In the opposite situation, the most efficient heating is the 
"DC" heating by direct current. In this case, the quasi-static 
energy input allows the accumulation of magnetic energy, the 
generation of currents and finally heating, for example, by 
reconnection of magnetic field lines dPriest & Forbes 20001 ) 
or by resistive dissipation of current sheets. However, in this 
simple version, these mechanisms are not efficient enough 
to explain coronal heating. Additional processes are thus 
generally introduced like turbulence, which can generate 
small scales where dissipation is much more efficient (see 
e.g. Gomez & Ferro Fontan 1988 Gomez & Ferro Fontan 1992 
Heyvaerts & Priest (1992") Einaudi et al. 1996 Dmitruk et al. 19971 



Galtier & Pouquet 1998| [Galtier and Bhattacharjee 2003| and 
Buchlin et al. 20071) . Moreover, turbulence could explain the 
measurements of nonthermal velocities revealed by the width 
of EUV and FUV lines. Observations of the transition region 
and corona of the quiet Sun from SUMER onboard SoHO 
reveal nonthermal velocities of about 30km s for temperatures 
around 3 x 10 5 K with a peak up to 55km s for some S IV lines 
dWarren et al. 1997[|Chae et al. 19981 ). 

Many observations of the solar atmosphere tend to show 
plasma in a turbulent state with a Reynolds number evaluated at 
about 10 12 . In particular, the most recent Hinode pictures seem 
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to show a magnetic field controlled by plasma turbulence at all 
scales (Nature 2007; see, also, IDoschek et al. 20071 1. Thus, the 
turbulent activity of the corona is one of the key issues to under- 
stand the heating processes. In the framework of turbulence, the 
energy supplied by the photospheric motion and transported by 
Alfven waves through the corona is transferred towards smaller 
and smaller scales by nonlinear coupling between modes (the so- 
called energy cascade) until dissipative scales are reached from 
which the energy is converted into heating. The main coronal 
structures considered in such a scenario are the magnetic loops 
which cover the solar surface, in active and quiet regions. Each 
loop is basically an anisotropic bipolar structure anchored in 
the photosphere. It forms a tube - or an arcade - of magnetic 
fields in which the dense and hot matter is confined. Because a 
strong guiding magnetic field (Bo) is present, the nonlinear cas- 
cade that occurs is strongly anisotropic with small scales mainly 
developed in the Bo transverse planes. In Figure Q] we present a 
schematic view of the turbulent energy spectrum expected for the 
coronal plasma which follows a power law over several decades. 
The larger scales may be determined directly by observation and 
correspond to about 10 5 km. The inertial range - that deter- 
mines the range of scales where the turbulent cascade operates - 
roughly starts between 10 4 km and 10 3 km, and extends down to 
unresolved (by spacecraft observations) scales which may be es- 
timated, from dimensional analysis, as of a few meters. The spa- 
tial resolution of different instruments is also reported to show 
the gap that we need to fill in order to completely resolve the 
heating processes that occur at the smallest dissipative scales. 
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Fig. 1. Schematic view of the energy spectrum of the coronal 
plasma in logarithmic coordinates. The spatial resolution of dif- 
ferent spacecraft instruments is reported: EIT/SoHO (1850 km), 
TRACE (350 km), SOT/Hinode (150 km) and Solar Orbiter (35 
km). The gray zone indicates currently unresolved small-scales 
which will be modeled by turbulent viscosities. 



The aim of the present study is coronal heating to perform 
of modeling by magnetohydrodynamic (MHD) turbulence. We 
explain the connection between the large scales (at which en- 
ergy is injected at loop footpoints through photospheric mo- 
tion) and the smallest scales (at which energy is dissipated and 
converted into heat). The input energy propagates through the 
network of magnetic field lines by "inward" and "outward" 
Alfven waves which nonlinearly interact and produce a turbu- 
lent cascade towards small scales. The foundation of our model 



is the one originally proposed by Heyvaerts & Priest (1992) 
where the unresolved small-scales are modeled through turbu- 
lent viscosities. These viscosities were extracted from an ad 
hoc EDQNM (Eddy-Damped Quasi-Normal Markovian) clo- 
sure model of MHD turbulence developed in spectral space for 
isotropic flows ( Pouquet et al. 19761. Spacecraft missions like 



SoHO, TRACE, or more recently Hinode show that this as- 
sumption is clearly not adapted to magnetic loops which are 
characterized by a strong longitudinal mean field Bo, compared 
to the magnetic perpendicular components, whose anisotropic 
nonlinear effects on the turbulent plasmas dynamics are more 
than likely, as many numerical simulations have shown. For that 
reason, in the present model, we use turbulent viscosities com- 
puted from an asymptotic (exact) closure model of MHD turbu- 
lence dGaltier et al. 20001lGaitier et al. 200"2l , also called Alfven 
wave turbulence. Signatures of such a regime have been detected 
in the middle magnetosphere of Jupiter (Saur et al., 2002). In 
our turbulent heating model, the unresolved small-scale equa- 
tions are perturbatively developed around a strong magnetic field 
Bo which leads to a strongly anisotropic turbulence. From the 
derived equations, we compute the turbulent viscosities which 
eventually allow us to obtain a self-consistent free-parameter 
model of coronal heating from which we predict a heating rate 
and a turbulent velocity that favorably compare with observa- 
tions. Figure [2] summarizes the schematic algorithm followed to 
obtain a self-consistent heating model. 
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Fig. 2. Schematic view of the self-consistent heating model. The 
resolved large-scales coronal flow is described by the "large- 
scale" MHD equations where turbulent (instead of standard) vis- 
cosity and resistivity are used to determine the inertial small- 
scale nonlinear dynamics. From these equations we compute the 
energy flux function released by the photospheric motion. The 
energy is transported by Alfven waves through the corona and 
the nonlinear interactions between wave packets generate an in- 
ertial cascade, ie. a production of smaller and smaller scales, that 
finally reach the dissipative smallest scales from which the en- 
ergy is converted into ohmic dissipation. The unresolved inertial 
small-scales are described by the asymptotic equations of Alfven 
wave turbulence from which the turbulent transport coefficients 
are computed. A self-consistent model is finally obtained by in- 
troducing the expression of these turbulent viscosities into the 
energy flux function. Then, the energy flux function is entirely 
determined and the heating rate can be evaluated. 

The organization of the paper is as follows: the next Section 
is dedicated to the large-scale description of magnetic loops 
where, in particular, we rederive some main results obtained by 
Heyvaerts & Priest (1992), using a more adapted notation. Then, 
in Section [3] the small-scale description is given and the turbu- 
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independent of the coronal plasma. More precisely, we assume 
for the boundary velocity field parallel to the y-axis: 



Fig. 3. Arcade geometry: the coronal plasma is confined in a vol- 
ume delimited by —h < x < +h, -oo < y < +00 and -I <z< +(■ 
The photosphere appears as two boundary planes at altitudes 
Z — -I and z = +€. Note the presence of Alfven waves prop- 
agating in opposite directions along Bq. 



lent viscosities are obtained. The model predictions (heating rate 
and turbulent velocity) for magnetic loops in active regions and 
in the quiet corona are given in Section |4] They are generalized 
to open magnetic lines in coronal holes in Section The last 
Section is devoted to the discussion and conclusion. 



2. Large-scale description of magnetic loops 

2. 1 . Geometry and boundary conditions 

We consider a set of magnetic loops - an arcade - anchored in 
the photosphere and subjected to the erratic motion of convective 
cells. For simplicity, this arcade has a rectangular cross surface 
which defines the perpendicular directions and it is elongated 
along the longitudinal magnetic field Bo = B(,e z that defines the 
parallel direction as shown in Figure [3] This magnetic structure 
is filled by a plasma which will be modeled as an incompress- 
ible MHD fluid. In the first part of this paper, we are concerned 
with magnetic loops; later on we will extend our results to open 
magnetic field lines for which the boundary conditions will be 
modified. For a magnetic arcade, the boundary limits are made 
of two planes at altitudes z-—( and z-+C which are connected 
by the uniform magnetic field Bo. Thus, the length of the loop in 
the longitudinal direction is 2€. The cross surface is delimited as 
—h < x < +h and -00 < y < +00. In other words, we assume 
an arcade characterized by a thickness of 2h, where h is the co- 
herence length of the magnetic field at the base of the corona, 
and with a large depth compare to the other length scales, with a 
translational invariance assumed in the y-direction. Note that we 
do not treat the expansion of the magnetic field at the footpoint 
loops where Bq has to be seen as an average value of the field at 
the base of the corona rather than in the photosphere itself. 

The energy reservoir for the solar corona resides in the sub- 
photospheric convective layer. Indeed, the high-/? plasma in the 
photosphere constrains the magnetic field lines to be convected 
by the plasma flow. These motions eventually induce a shear- 
ing and a torsion of coronal magnetic field lines. The boundary 
motions at footpoint loops are assumed to be perpendicular and 



u(x,z = €)= -u(jc,z = -€) = U(x)e 



> ■ 



(1) 



This input of energy is then propagated through the coronal 
medium where the low-/? confines the plasma into the magnetic 
arcade. These boundary motions are then able to produce and 
maintain a state of full turbulence in the corona. Note finally the 
absence of mass flow through the boundaries. 

2.2. The MHD model 

The standard incompressible MHD equations are used as a first 
approximation to model the solar plasma, namely: 



d,u + u Vu 

V u 
<9,B 

V B 



-VP/po + j x B/po + vAu, 
0, 

V x (u x B) + 77AB, 

0, 



'"7 + £j + 



S-D, 



w -\ E x B 

T = u-\ Po - + P-cr\ 



E + u x B = p rjj, 



(2) 
(3) 
(4) 
(5) 

(6) 

(7) 
(8) 



where u is the velocity, B the magnetic field (such as B = Bo+b), 
j = V x b/po the current density (where po is the permeabil- 
ity of vacuum), E the electric field, & the viscous stress tensor 
(with o-jj = pQv{dui/dr j + duj/dr,)), v the kinematic viscosity, 
r\ the magnetic diffusivity and po the uniform mass density, and 
P is the gas pressure. In the energy equation ©, the right hand 
side term S is the large-scale energy input created by the photo- 
spheric motion and D is the dissipation which mainly happens 
at the smallest scales because of, e.g., viscous and resistive ef- 
fects. The left hand side term, V • T, denotes the nonlinear trans- 
fer of total (kinetic plus magnetic) energy. (Note that all terms 
in equation © are expressed in unit volume.) Equation (0 de- 
scribes how the total energy is nonlinearly transferred at differ- 
ent scales by velocity transport, by work due to the gas pres- 
sure, the viscous stress or the electromagnetic work due to the 
Poynting vector (E x B/po). Other (higher order) processes like 
thermal conductivity which happens at microscopic scales are 
not considered here. Finally, equation © is the standard Ohm's 
law which allows us to close the system. 

The basic idea of this coronal model is the separation of 
the non-resolved inertial small-scale dynamics (see FigureQ]and 
lU from the large-scale dynamics detected with current instru- 
ments. In practice, we assume that the previous MHD equations 
describe the large-scale behavior of the plasma; the nonlinear 
small-scale processes (ie. the effects of the inertial small-scales 
on the large-scale dynamics) will be modeled through effective 
dissipations, also called turbulent viscosity and resistivity. In 
other words, the usual molecular viscosity and magnetic resistiv- 
ity are replaced by turbulent coefficients which are several orders 
of magnitude larger. The expression of these turbulent viscosities 
(hereafter denoted v, and rj,) will be derived in Section[3]and the 
MHD equations in the next Section. 

An important assumption to derive this heating model is that 
the slow boundary motions are able to sustain a fully turbulent 
state in the corona, with a typical photospheric time longer than 
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the coronal time scales. Basically, this statement means that the 
characteristic time of the photosphere has to be longer than the 
time needed for a perturbation to cross the coronal loops. For 
a typical longitudinal direction I of 10 7 m and an Alfven speed 
of several 10 6 m/s, one finds a crossing time of the order of Is. 
For a photospheric granule with a size A = 10 6 m and a veloc- 
ity u = 10 3 m/s, one finds a lifetime of about 10 3 s. In this case, 
the time scale separation is thus well satisfied. However, a valid- 
ity condition may arise if the turbulent nature of the photospheric 
motion is taken into account since their velocities seem to follow 
a Kolmogorov spectrum dRoudier et al. 19871 IChou et al. 19911 



Espagnet et al. 1993 1. Then the photospheric velocity at a given 
scale A should scale as A 1 ^ 3 and the typical photospheric nonlin- 
ear time as A 2 ' 3 , which decreases with the size of the granules. 
The associated condition is discussed and used in SectionP 



2.3. Large-scale solutions 

The nonlinear solutions of the previous MHD equations are 
non trivial and require direct numerical simulations on super- 
computers. However, we are here mainly interested in the large- 
scale behavior of the coronal flow which allows us to make some 
simplifications. We will seek solutions of the form 



U = u(x,Z)£y, 

B = B e z + b(x,z)e y , 



(9) 
(10) 



independent of y and f, which satisfy, in particular, the diver- 
gence free conditions (01 and These solutions do not depend 
on time since we assume magnetic loops maintained in a station- 
ary regime of turbulence. We also observe that the fluctuating 
velocity and magnetic fields are taken perpendicular to the uni- 
form magnetic field Bo. This assumption is compatible with the 
fact that, in strongly magnetized media, pseudo-Alfven waves 
are mainly along the uniform magnetic field whereas the nonlin- 
ear dynamics is dominated by shear- Alfven waves in the trans- 
verse planes. We recall that shear-Alfven and pseudo-Alfven 
waves are the two kinds of linear perturbations about the equi- 
librium, the latter being the incompressible limit of slow mag- 
netosonic waves. Replacing the previous solution types into the 
MHD equations © and © leads to 

d x (P + b 2 (x,z)/2» Q ) = 0, (11) 

d z (P + b 2 (x,z)/2 Mo ) = 0, (12) 

Bod z b(x,z) + [io pov,Au(x,z) = 0, (13) 

B d-u(x,z) + T],Ab(x,z) = 0. (14) 

These equations describe the large-scale evolution of the plasma 
once the turbulent regime is established. In particular, they al- 
low us to calculate the pressure field when the magnetic field 
is known. The linear form of equations ( fT3l and ( TBI may be 
misleading since the small-scale nonlinearities are concentrated 
in the turbulent coefficient transports, namely the turbulent eddy 
viscosity v, and magnetic diffusivity tj,. This assumption can be 
justified when looking at the nonlinear evolution equations for 
the energy spectra, as used in Section[3] Indeed, in Fourier space, 
a nonlocal analysis leads to diffusive effects due to the small- 
scale dynamics over the large scale motions, actually stemming 
from absorption terms involved in the modeled energy transfers 
between small and large scales. 

2.4. Photospheric conditions 

The photospheric forcing introduced in equation (Q]i mimics the 
large-scale motion of magnetic footpoints imposed by the con- 



vective cells. Following |Heyvaerts & Priest (1992.)} we assume 
that the frozen-in law is satisfied in the dense photosphere which 
leads, when coupled with the Ohm's law ©, to the boundary 
relation for the electric field: 



E(x,z = ±€) = +U(x)B e K . 



(15) 



In the close neighborhood of the photosphere, this relation is 
modified to take into account the magnetic diffusivity. This gives 
for the x-component: 



db 

E x (x,z - +£) = +r] t —(x,z) + U{x)B(). 
dz 



(16) 



The continuity of the tangential component of the electric field 
implies that 



db 

—(x,z = ±e) = 0. 
dz 



(17) 



These photospheric conditions are of course compatible with 
equations ( fT3l and ( TBI when z - +1. 



2.5. Energy flux injection 

The aim of this section is to express the energy flux released by 
photospheric motion along coronal loops. We are thus interested 
in the z-component of the flux (0 at the boundaries z = ±(, 
namely: 



r z (x,z = ±€)=e z 



u _\ ExB 

po— +P-& + 

2 / Ho 



(18) 



The velocity boundary conditions (HJ do not involve (convective) 
mass flows through the photospheric planes since u z (z = +1) = 0. 
Therefore, the two first terms in expression ( fT8l do not con- 
tribute, and the input flux reduces to 



T z (x,z = +1) = Ti z +r z , 



with 



n z = - ( u • <x)_ , p z = 



ExB 

Mo 



(19) 



(20) 



i.e. only the viscous stress tensor and the Poynting vector con- 
tribute to the energy flux transfer to the corona. One of the main 
goal of this paper is to evaluate this energy flux by first solving 
equations ( fT3l l and (fT4l i. which depend on the turbulent viscosity 
and resistivity, and second, by calculating these turbulent coeffi- 
cients from Alfven wave kinetic equations. 

2.6. Solving the large-scale dynamics 

In this section, we solve equations $13[ and ( fT4b to find the large- 
scale plasma behavior. We develop the x-dependence of the ve- 
locity and magnetic fields in Fourier series (since h is assumed 
to be the coherence length of the magnetic and velocity fields in 
the x-direction) 



£\z) cos — + 2, £f(z) sin — , 



(21) 



where £ is either u or b. Note that other solutions are pos- 
sible with, eg. a twisted flux tube ( jlnverarity and Priest 1995[ >, 
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which do not drastically change the coronal heating predictions. 
Substitution of ( 1211 into < fl~3T > and ( [T4l gives 



= 0, 



dbf I d 2 n 2 n 2 

duf Id 1 n¥\ lffl 
where (;) stands for indices (1) or (2), and finally 



(22) 
(23) 



4„(0 



d u 



dz 4 



I^ +2 (^)Y^ + ('f) 4 u^0. (24) 
HoPoViT], \ h I J dz L \ h I 



Note that the magnetic field satisfies the same equation, only the 
boundary conditions make the distinction between the two fields. 
Using the following notation: 



B aA = rm iJn _ d m 
^JxopWtVt ' " h ' d(az)" 

equation (l24l rewrites 

r 4 M »-(i+2^„V"j; ) +4"!; ) = o, 



(25) 



(26) 



whose solutions may be expressed in terms of exponentials. The 
photospheric boundary conditions (Ql and (TTTl i combined with 
equation (l22l gives the relations 



uf(z 



,1 (i) 

d-u,/ 



-.U® and ^(z = ±0 = ±a 2 ^U®, 
dz- 



which lead eventually to the solutions 

u%) = 



%- sinh(3t>z) %+ srah^az) 



sinh(%+a£) 



where 



<x> ± _ 



VI +4^+1 



(27) 



(28) 



(29) 



The substitution of (f28b into d22l and d23l gives, after some 
manipulations, the following expression for the magnetic field 
Fourier coefficients 



bf(z) = 



BoU),' } (%+ cosh(%- az) R~ cosh(2l>z) 



K + + %~ \ rja smh(%-ai) ?ja sinh(%+at 



(30) 



2. 7. Evaluation of the energy flux 

From the solutions d28l and (130b it is now possible to calculate 
the energy flux ( fT9b at footpoint levels z = +£. The contribution 
of the Poynting vector to the flux is 



P z (x,±€) = + — U(x)b(x,z = ±t), 
Ho 

and the contribution of the viscous stress tensor is 

du(x, z = +() 



U z (x,±{) = +p v,U(x)cr~ y = +p v,U(x)- 



dz 



(31) 



(32) 



The input flux coming from the photospheric boundaries (z = 
+€) is calculated by averaging over the x-periodicity of the ve- 
locity and magnetic fields ( 1211 : 



\f z (z = ±01 = ^ 

2/2 



[P z (x, +{) + n z (x, ±€)] dx 

h 



(33) 



After some manipulations, we obtain the total flux (from the two 
footpoint levels): 



\r z \ = p £ 



B ° -u 2 

- w n 



sinh (q;^ 



k cosh [ +4^a^J - cosh (at 



(1 + 2^) sinh [ VI + 4Alat] / y/l +AA 2 n " 
cosh 

[Vl + 4Alat] - cosh (al) 



(34) 



with Ul = C/ (1 >„ + J/ (2) „. Note the appearance of a turbulent 
magnetic Prandtl number (v,/tj,) whose origin is directly linked 
to the action of nonlinear small-scale processes on large-scales. 
We will later show that this number may be taken equal to unity 
in the Alfven wave turbulence regime. Also note that the turbu- 
lent viscosities enter in expression ( f34b through the coefficients 
a and A„. The origin of the presence of v, and rj, is thus due, on 
the one hand, to the large-scale MHD equations (fT~3b — (fT~4-b and, 
on the other hand, to the Poynting vector OTT l and the viscous 
stress tensor (1321 which also include the small-scale nonlinear 
retroaction. The energy injected from the magnetic footpoints is 
propagated through the corona by Alfven waves. Then, the col- 
lisions between "upward" and "downward" wave packets lead 
to nonlinear dynamics which is characterized by a turbulent cas- 
cade that transfers energy from the energy injection scales down 
to the smallest dissipative scales where energy is destroyed, and 
the corona heated. The coronal loops reach a stationary state of 
fully developed turbulence when the dissipation equals the elec- 
tromagnetic and viscous stress. 

2.8. Turbulent spectrum of photospheric velocities 

The coronal heating flux is now evaluated from d34l . The tur- 
bulent viscosities entering into this flux expression will be 
computed in the next Section. Therefore, the last unknown is 
the magnitude of the velocity at the photospheric boundaries 
(namely U„). The goal of this section is to evaluate the scale 
dependence of photospheric velocities. 

As reported by different measurements ( IRou dier et al. 19871 
IChou et al. 19911 Espagnet et al. 1993 1, the photospheric veloc- 
ities are assumed to follow a Kolmogorov isotropic spectrum 



22,2/3,-5/3 



with 



u 2 = I u 2 (k)dk = u ll 2 



(35) 



(36) 



where u is the photospheric rms velocity, U 2 is the power spec- 
trum of the boundary velocity fluctuations, and ki„j is the largest 
scale at which energy is injected into the system; it corresponds 
to the size of the largest granular cells, i.e. about 1000km. In 
other words, it is the integral scale of the turbulent photospheric 
flow from which the inertial range starts. The wave number 
kj n j may be directly connected to the thickness of the arcade, 
such that ki„j = njh. More generally, the wave number k and 
the harmonic number n may be connected through the relation 
k = nnjh. Converting the sum in equation ( 1361 into an integral 
gives the relation 



Ul = \u 2 (k) = \u 2 nV\ 
h 3 



(37) 



The basic idea of the model is that Alfven wave packets in- 
teract nonlinearly along coronal magnetic structures in which the 
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plasma is assumed to be maintained in a fully turbulent state be- 
cause of the slow photospheric motion. We have seen in Section 
12.21 that a condition may arise on time-scales since the photo- 
sphere is also turbulent and its typical nonlinear time, r^, may 
be very short for small granules. In order to satisfy this condi- 
tion, we impose that t'^(A) > r cross , where T cross is the time 
needed for a perturbation to cross half of the loops of length It 
at speed Bo I -\fpojJo, namely 



Bq ' 



and t^(A) is the eddy turnover (nonlinear) time of granular 
structures at a typical scale A and velocity u(A); 

A A 2 ' 3 /2n\ l/3 



< h M) 



u(A) u \k inj 



(39) 



This time-scale separation induces for the wave number k 
2n/A, the inequality k < k cross , where 

.3/2 



7T I h B / VPo/"o\ 

kcross = tN crass, With N cross = 12- - I . (40) 

Hence, the upper limit of the summation on harmonic numbers 
in d34l i is given by N cr oss which is a large number since the 
Alfven velocity is much larger than the photospheric rms ve- 
locity. For that reason, the expression (l34l may also be written 
as an integral with 



(41) 



The energy injection produced by the photospheric motion is fi- 
nally rewritten as follows: 



\r z \ =p 

with 
¥(*) = 



/vT_Bo 



sinh(o^ 



U 2 {ky¥{k)- dk, 



(42) 



cosh [ Vff 2 + 4k 2 t] - cosh(atf) 

a 2+2k 2 sinh[Va 2 +4fc 2 l] 
a Va 2 + 4k 2 cosh [ Va 2 + 4A: 2 ^] - cosh(crf) 



(43) 



Under the time-scale condition discussed above, the coronal 
structure is submitted to slow external forcing. Note that other 
time scales are present in our problem; namely, the nonlinear 
coronal time, t„/, and the Alfven time, ta, which has nothing 
to do with the crossing time r crosj introduced above. It may 
be interpreted as the time of interaction between two counter- 
propagating Alfven wave packets. Both times will be used in 
Section[3]where the nonlinear small-scale dynamics will be an- 
alyzed. The main assumption made in such wave turbulence dy- 
namics is that the Alfven time is smaller than the nonlinear coro- 
nal time. Note that in MHD the transfer time-scale is built on t„; 
and ta which makes a great difference in Navier-Stokes fluids 
for which the transfer time is the nonlinear time. This point will 
be further discussed in the next section. Then the following hier- 
archy of times may be estimated 



As seen in Section 12.21 the nonlinear photospheric time can 
be as large as 10 3 s and the crossing time is about Is. For 
the two other characteristic time scales, no observational con- 
straints can be used. However, recent direct numerical simu- 



lations (Bigot et al. 2008b i show clearly that a turbulent MHD 



flow evolving in a strong uniform magnetic field (Z?o > b) ex- 
hibits a time ratio ta I t,,i < 1 > which decreases as the strength of 
the uniform field increases. 



(38) 3. Small-scale dynamics and anisotropy 



In the previous Section, we have been able to evaluate the flux 
of energy released by photospheric motion along coronal loops. 
The calculation has been greatly simplified by the assumption 
of a scale separation between large and small (inertial) scales, 
and thus the introduction of turbulent eddy diffusivities (v, and 
T],). In fact, all the complexity of the nonlinear dynamics is con- 
centrated in these coefficients. We will now try to evaluate them 
by going back to previous works on Alfven wave turbulence. We 
first summarise previous results on anisotropic MHD turbulence. 
We then introduce the asymptotic equations of Alfven wave tur- 
bulence which are not rederived here (see [Galtier et al. 20001 >, 
and show their main properties. Finally, we derive the turbulent 
viscosities from such nonlinear equations. 

3.1. Role of anisotropy in MHD turbulence 

The first description of incompressible MHD turbulence 
proposed by Iroshnikov and Kraichnan (IK) in the 1960s 
(Iroshnikov 1964, Kraichnan 1965) is based on a modified ver- 
sion of the Kolmogorov dimensional analysis where the large- 
scale magnetic effects are taken into account although the fluid 
is supposed to be globally isotropic. The IK's heuristic model 
includes Alfven waves which approximate the magnetic field at 
the largest scales as a uniform magnetic field. Then, it is the 
sporadic and successive collisions between counterpropagating 
Alfven waves that produce the turbulent cascade. At leading 
order, the nonlinear dynamics is driven by three wave interac- 
tions, and the transfer time through the scales is estimated by 
Ttr ~ t Ii/ t a, where t„\ ~ A/ux is the nonlinear turnover time 
whereas ta ~ A/Bq is the Alfven wave period (where Bq stands 
for the Alfven speed, taken as the rms magnetic field fluctua- 
tion). Hence, the energy spectrum is k~ 3 ^ 2 . 

The weakness of the IK phenomenology is the apparent con- 
tradiction between the presence of Alfven waves and the ab- 
sence of a strong uniform magnetic field. This point and the 
role of a strong external magnetic field Bo have been widely dis- 
cussed in the literature, in particular during the last two decades 
(Strauss 1976, Montgomery & Turner 1981; Shebalin et al. 
1983; Oughton et al. 1994; Goldreich & Sridhar 1995; Ng & 
Bhattacharjee 1996; Matthaeus et al. 1998; Galtier et al. 2000; 
Nazarenko et al., 2001; Milano et al. 2001; Chandran 2005; 
Bigot et al. 2008a). One of the most clearly established results 
is the bi-dimensionalization of a turbulent flow with a strong re- 
duction of the nonlinear transfers along the Bo ambiant field. 

In the case of coronal structures, like magnetic arcades 
or open magnetic fields in coronal holes, a strong magnetic 
field at large-scale is clearly present. This field is often mod- 
eled as a uniform field through, for example, the approxima- 
tion of reduced MHD (see, e.g., Dmitruk & Matthaeus 2003; 
Buchlin et al. 2007). It was thought until recently (see eg. 



TA < T„l < Tcross < T 



(44) 



Oughton et al. 2004 1 that the reduced MHD approximation was 
only able to describe the regime of strong turbulence but a re- 
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cent work (Nazare nko 20071 ) has extended its domain of appli- 
cation to the wave turbulence regime as well. This point is par- 
ticularly important for our model which is based on the latter 
regime. In our model, the magnetic arcade is composed of a net- 
work of parallel strands that are maintained in a fully developed 
turbulent regime. Contrary to the isotropy hypothesis originally 
assumed by Heyvaerts & Priest (1992) we include anisotropic 
effects through turbulent viscosities. In this approach, it is thus 
necessary to distinguish between the perpendicular (±) and par- 
allel (||) directions to the uniform Bo field. 

An important anisotropic property discussed in the literature 
(Higdon 1984; Goldreich & Sridhar 1995; Galtier et al. 2005; 
Boldyrev 2006) is the interdependence of perpendicular and par- 
allel scales. According to direct numerical simulations (see, e.g., 
Cho & Vishniac 2000; Maron & Goldreich 2001 ; Shaikh & Zank 
2007; Bigot et al. 2008b), one of the most fundamental results 
seems to be the (critical) balance between the nonlinear and 
Alfven times which leads to the scaling relation k\\ ~ k 2 / 3 . In 
other words, a turbulent MHD flow evolving in a strongly mag- 
netized medium seems to be characterized by an approximately 
constant ratio (generally smaller than one) between the Alfven 
and the nonlinear times. This result (generally) implies a dynam- 
ics mainly driven by Alfven wave interactions. 

3.2. Alfven wave turbulence 

The incompressible MHD equations (f2])-(|5]l may be reformu- 
lated in terms of the Elsasser fluctuations 



u + sb . 



(45) 



with magnetic fields written in velocity units (b = b/ -\fjJoPo), 
and in the inviscid case (ie. in the absence of the dissipative terms 
proportional to the kinematic viscosity and the magnetic diffu- 
sivity), 



d,z s - sB ■ Vz s 
V ■ x s 



-z 
0. 



Vz s - VP* , 



(46) 
(47) 



Here, s — +, stands for the directional polarity, indicating the 
propagating direction of waves along Bo = Bo/ y/zopo; P* is the 
total (kinetic plus magnetic) pressure. Note that the present anal- 
ysis is focused on the nonlinear plasma dynamics from which 
the turbulent transport coefficients will be derived. In the pres- 
ence of a strong uniform magnetic field, MHD turbulence may 
be dominated by wave dynamics for which the nonlinearities are 
weak, and thus ta <k t„/. In this limit, a small formal parameter 
e can be introduced to measure the strength of the nonlinearities 
to give, for the jth-component 



(d, - sB d\\)Zj = -ez m s d m z S j - djP m 



(48) 



(with Einstein's notation used for indices). Note that the parallel 
direction (||) corresponds to the z-direction. 

We will Fourier transform such equations, with the following 
definitions for the Fourier transform of the Elsasser field compo- 
nents z"j(x, t): 



(k, t) e 



c/k, 



(49) 



where <±>k — Bo^y is the Alfven frequency. The quantity a*(k, f) is 
the wave amplitude in the interaction representation, hence the 
factor e' SWlt . The Fourier transform of equation d48l ) thus is 



cu;.(k) 



-i€k m P j n 



I 



a- m s {K)a s n {lS)e isAu,t 8^dKdL 



(50) 



Here, Pj n is the projector on solenoidal vectors such that 
P/„(k) = 6j„ - kjk n /k 2 ; (J^l = #(k - k — L) reflects the triadic 
interaction, and Aa> = <±>i - u)k - <i) K is the frequency mixing. 
The appearance of an integration over wave vectors k and L is 
directly linked to the quadratic nonlinearity of equation d48l (as 
a result of the Fourier transform of a correlation product). 

Equation (T50b is the compact expression of the incompress- 
ible MHD equations when a strong uniform magnetic is present. 
It is the point of departure of the wave turbulence formalism 
which consists of writing equations for the long time behav- 
ior of second order moments. In such an asymptotic and sta- 
tistical development, the time-scale separation, ta/t„/ «: 1, 
leads to the destruction of some nonlinear terms, including the 
fourth order cumulants, and only the resonance terms survive 
dGaltier et al. 20001 iGaltier et al. 20021 . It allows one to obtain a 
natural asymptotic closure for the moment equations. In such a 
statistical development, the following general definition for the 
total (shear- plus pseudo- Alfven wave) energy spectrum is used; 



<a*(k)a*(k')> = E s (k) 6(k + k')/k x . 



(51) 



where () stands for ensemble average and k ± = |kj. In absence 
of helicities and in the case of an axially symmetric turbulence, 
the asymptotic equations simplify. For the shear-Alfven waves, 
the energy spectrum is given by 



(52) 



where f(k\\) is a function fixed by the initial conditions (ie. there 
is no energy transfer along the parallel direction). The transverse 
part obeys the following nonlinear equation (the small parameter 
e is now included in the time variable) 



d,E s JkJ = J- J Jcc 

[k ± E s ± (L ± ) - L ± E[(k ± )] d K± dL ± , 



k± 



(53) 



where <p is the angle between k ± and L ± , and 9 is the angle be- 
tween k ± and k± with the perpendicular wave vectors satisfying 
the triangular relation k ± = h± + k ± (see Figure |4]i. Note that 
from the axisymmetric assumption, the azimuthal angle integra- 
tion has already been performed, and we are only left with an 
integration over the absolute values of the two wave numbers, 
k ± = | if j. | and L_l = |LJ. In the same way, equations can be 
written for pseudo- Alfven waves which are passively advected 
by shear-Alfven waves, namely 



d t El(k ± )=^- f f sin0^£7> ± ) 

[ifcjEjfCLj.) - L x E^(k x )] d K± dL ± , 
with by definition 

Kseudo(k±A)=m)EZ(k±), 



(54) 



(55) 



where f(k\{) is a function fixed by the initial condition. Equations 
(|53l and (1541 - called wave kinetic equations - only involve 
perpendicular nonlinear dynamics (with a dependence only on 
the perpendicular components of the wave vectors), a situa- 
tion expected from previous works (Strauss 1976, Montgomery 
& Turner 1981; Shebalin et al. 1983; Oughton et al. 1994; 
Goldreich & Sridhar 1995; Ng & Bhattacharjee 1996; Matthaeus 
et al. 1998) that showed the decrease of the parallel nonlinear 
transfers with the increasing strength of the external magnetic 
field B . 
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BJ/z 




Fig. 4. Triadic interaction k = k + L and its projection in the 
plane perpendicular to Bq (such that k ± = k± + L ± ). 



3.3. Anisotropic turbulent viscosities 

The main goal of this Section is the derivation of turbulent 
viscosities in the context of strongly anisotropic MHD turbu- 
lence. We have seen that wave turbulence leads to the inhi- 
bition of nonlinear transfers in the Bo direction, and to the 
bidimensional reduction of the dynamics. This result makes 
a strong difference to purely isotropic MHD turbulence from 
which the turbulent viscosities were derived in the original 



model (Heyvaerts & Priest (1992)). We believe that taking into 



account the flow anisotropy is an important improvement in the 
description of the coronal heating problem. 

The wave kinetic equations are assumed to describe the non- 
linear dynamics at small inertial scales. These scales correspond 
to the scales unresolved by the current spacecrafts which are re- 
ported in Figure Q] The turbulent viscosities represent the aver- 
age effect of nonlinearities at these small scales, and they appear 
in the large-scale equations (fT~3T> — (TT~4T>. In practice, we make a 
Taylor expansion of the wave kinetic equations for nonlocal tri- 
adic interactions (Pouquet et al. 1976 1 that satisfy the relation 



which leads to 



L ± - k± - k ± cos 9, 



(56) 



(57) 



for triadic interactions k ± = k x + L ± (see Figure [?}. Then, re- 
lation d57l ) is substituted into the kinetic equations d5"3l-(l5"4li for, 
respectively, the shear- Alfven waves 



J DO JO 



d t E s x (k ± ) = -klEi(kJ 

In 



cos 2 sin 2 6d6, 



(58) 



and the pseudo- Alfven waves 

d t E*(k ± ) = -J£j5jf(*j.) 

f ^-El s (K ± )dK ± f sin 2 OdO. 
J Bq Jo 

Finally, we obtain: 

d,E[(k ± ) = -v{kiE s x (k ± ), 
d t E°,(k ± ) = -vjtfEJOkJ, 



(59) 



(60) 
(61) 



where and vjj are the so-called turbulent viscosities, corre- 
sponding respectively to shear- and pseudo-Alfven waves, de- 
fined as 



2 r- 

v\ = ^- El\K ± )dK x . 
BO J 



(62) 
(63) 



The turbulent viscosities are thus directly obtained after inte- 
gration over the 6 angle. Note that, similarly to the original 
paper (Heyvaerts & Priest 1992), the integration is made from 
the beginning of the inertial range (at Ki„j) where the injec- 
tion of energy is made; this choice contrasts with some models 
where a wavenumber cut-off is chosen inside the inertial range 
which gives a k-dependence of the turbulence viscosity (see eg. 
Baerenzung et al. 2008). An evaluation of the turbulent viscosi- 
ties is now possible by the use of the exact power law solutions 
of the integro-differential equations d53l and d54l ). This is first 
made in the context of a magnetic arcade for which the inward 
and outward Alfven waves are exactly balanced (balanced case). 
This approach is then generalized to open magnetic fields for 
which outward Alfven waves are dominant (unbalanced case). 

3.4. Balanced turbulence for a magnetic arcade 

For balanced turbulence, as said before, it is not necessary to 
distinguish outward from inward Alfven waves. We thus drop 
the directional polarity s, and we obtain 



2 /- 

tt 2 r 



_(K±)dK ± . 



(64) 
(65) 



In this case, the exact power law solution of equation d53l l is 
dGalfier et al. 2000b 



(66) 



where Ck is the Kolmogorov constant whose value is about 
0.585, and P ± is the constant energy flux of shear-Alfven waves 
whose definition is 



d t E ± (k x ) = -d k± P ± (kJ . 



(67) 



It is important to note that it is the condition of constant en- 
ergy flux that allows to find (by the Zakharov transformation) 
the exact power law solution of the nonlinear wave turbulence 
equation. The introduction of (|66l l into d64l ) and d65l l gives 



(68) 




(69) 



Note the relation vy = 4v ± between the turbulent coefficients 
corresponding to pseudo- and shear-Alfven waves. In particular, 
this means that the pseudo-Alfven wave energy is susceptible to 
decay faster than the shear-Alfven one. Actually, this property 
has been detected recently in numerical simulations of freely de- 
caying MHD turbulence (Bigot et al. 2008a). Finally note that in 
this wave turbulence regime, there is an automatic equipartition 
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state between the kinetic and magnetic energies. Therefore, no 
distinction will be made between the average nonlinear effects 
acting on the velocity and on the magnetic field. A unit turbulent 
magnetic Prandtl number will then be taken (ie. v,lrj, = 1). 



4. Magnetic arcade heating 

4.1. Connection between large- and small-scale dynamics 

The magnetic arcade being supposed in a fully turbulent state, 
the energy transfer rate is then assumed to be constant at each 
scales in the inertial range, from the largest photospheric scales 
at which energy is injected, down to the small dissipative scale 
fom which heating starts to occur. This fundamental remark al- 
lows us to link the first estimate of the energy flux ( l42l i. based 
on the photospheric motions, to the flux through the anisotropic 
arcade scales, estimated from the turbulent eddy viscosity and 
magnetic diffusivity. 

According to relations d52l ) and d67| i the energy flux P ± , in- 
troduced above, satisfies the general relation 



d t E shear (k x ,k\\) = -f(k\\)d k± P x . 



(70) 



We remind that / is an arbitrary dimensionless function which 
reflects the absence of nonlinear transfers along the Bo direction. 
The coronal heating in a fraction of the arcade, say Q, can thus 
be estimated from the small-scale dynamics by 



Q = 2(2hLp P x 



j f(k)dk h 



(71) 



where L delimits a fraction of the arcade thickness in the y- 
direction and Q is measured in Js~' (Note the presence of the 
constant mass density). In the same manner, we may evaluate 
the coronal heating from the large-scale energy flux \T : \. It gives 
for the same coronal volume: 



Q = 2hL\T z \. 



(72) 



Hence, the relation between both photospheric and turbulent 
fluxes: 



P± 



(73) 



4.2. General expression of the heating flux 

The small-scale analysis gives two turbulent viscosities coming 
from, respectively, the shear- and pseudo-Alfven wave energy. 
To be consistent with our model, we will only consider the for- 
mer one since, in the strongly anisotropic limit of wave turbu- 
lence, the pseudo-Alfven waves fluctuate only in the parallel di- 
rection to the ambiant magnetic field. However, we note that the 
inclusion of the pseudo waves effect modifies only slightly the 
form of the turbulent viscosity (by a factor 5). Then, relation (TTBl 
gives (with v, = v x ) 



n 4 ' C\ 2£ P0 ffWh 



(74) 



Together with the injection wave number K in j - ir/h, related to 
the arcade thickness, we obtain 



32 



n 2 C\ h 2 



fa 



'oB J 



f(k\\)dk n . 



The substitution of ( l42l into d75l ) leads to the relation 

32 v 2 ( r rK, 



[ f{ky)dk\\= f"°"u 2 (kmk)dk, (76) 

J Jk,,,: 



where a unit turbulent magnetic Prandtl number is taken (v, = 
T] t ). The above expression also gives 



f*cnm 

A a 2 ^{k)k- 5l3 dk = 1, 

J km; 



where 



An = 



nC 2 /7T\ 5 / 3 h 3 u 2 1 



T\h) IE 



24 \h) (B\ f m)dh 
In the next section, we use the following approximation: 



(77) 



(78) 



(79) 



4.3. Different scale contributions 

According to the specific value of the k wave number, different 
estimates may be made for the function ^(k) that enters in the 
expression ( l42b of the energy injection by photospheric veloci- 
ties, namely 



a <tz k 



a 



« k « a : *¥ 2 (k) 2, 



k <sc 



*3(*)^. 



(80) 
(81) 
(82) 



We recall that a = Bo/ -^fiopoVii], = Bq/v,; it is therefore a way 
to measure, for example, the relative importance of the diffusive 
terms in equations (Tl3T > and (TBI) which can be rewritten as 



ad z b(x, z) + Au(x, z) = 0, 
ad z u(x, z) + Ab(x, z) = 0. 



(83) 
(84) 



On can also more easily evaluate the relative importance of the 
Poynting vector OTb and the stress tensor d32l to the flux with 
the relation 



f> z (x, +€) + n z (x, +€) = 
+poU(x)v t [ah(x,z = +() + d z u(x,z = +()]■ 



(85) 



For the case a <s k, it is clear that the main contribution to the 
flux comes from the stress tensor, whereas for the two other k in- 
equalities, the main contribution comes from the Poynting vec- 
tor. Relation ( TTTb can thus be written in these different regimes 
as: 

(0 a «: k inj , 



A ar' 



[ kC '°' S ¥ 1 (k)k- 5/3 dk= 1; 

•Jk,,,, 



(86) 



(it) k, 



mj « a « ki/ 



(75) A a 2 



' n a (Across 

I x V 2 (k)k- 5/3 dk+ I Vdtyk-^dk 



(87) 
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(Hi) k inj C <K a <K kcross, 



Aq a 2 



x ¥ 3 (k)k~ 5/3 dk + V 2 

k mj J VuTf 



(k)k- 5/3 dk 



j"'"\ l (k)k- 5/3 dk\ 



(88) 



0'v) k cross « a «: k 2 cro J, 



x ¥ i (k)k- 5/3 dk+ *¥ 2 {k)k- 5l3 dk 

k„„ J^IaJl 



= 1; 



(89) 



(V) KroJ « 



A cr 



(a,k)k- 5,3 dk = 1. 



These relations simplify respectively to: 
(/) 3Aoa(^ 3 «-^ 7 3 ) = l; 

(b) 3A (o*^„ + a\ 2/3 - 2a 4/3 ) = 1; 



(nO 3Ao ^ + ~^ 3 + " 2« 4/3 ) = 1; 



(90) 

(91) 
(92) 

(93) 



(fv) 3A f -/tr 8/3 - a 2 k-^ + Za 5 ' 3 * 1 ' 3 ) = 1 ; (94) 
l 8£ J 8 / 



M Iao^K /3 -^ 3 ) = l 



(95) 



These five different regimes correspond to a different intensity 
of the small-scale nonlinearities. This intensity is measured in 
terms of an equivalent wave number a and it is compared to the 
two available scales, namely kj„j and k cross . Case (i) is the situa- 
tion where the small-scale nonlinearities are the strongest: they 
are felt beyond k m j. Case (v) is the situation where the small- 
scale nonlinearities are the weakest: the diffusive terms in equa- 
tions d83l-(l84ll are then negligible. 



4.4. Predictions for a coronal plasma 

We define the following basic coronal quantities: 

h = 10 6 Wm, 
t = IQ 1 L m, 
Bo = l(T 2 So T, 



ft = 10 3 1/ ms" 1 , 

p = l(T 12 A4kgnr 3 , 

which give the estimates: 

B = 8.9 x 10 6 S AC 1/2 ms" 1 , 

k inj = 3.1 X lO^TT 1 nT 1 , 

k cross = 2.4xl0- 4 -H 1/2 £ 3/2 'Z/ 3/2 S 3/2 A1- 3/4 m- 1 , 

A = 1.2xl0"'W 4/3 1/ 2 S 2 Mm 4/3 . 



(96) 
(97) 
(98) 
(99) 
(100) 

(101) 
(102) 

(103) 

(104) 



With these values, the only satisfied equation is d94b which cor- 
responds to the inequality of case (iv). The solution of equation 
(IH is: 



a = 4.1 x 10^ 



-4/3 ,7,2/3 



AT 1/3 W 



-2/3 



(105) 



Then, it gives the following estimate for the turbulent viscosity: 
v, = 2.2 x 10 9, W 4/3 'Z/ 2/3 £ 73 AT 1/6 mV 1 . (106) 
Hence, the heating flux per unit area 

\T Z \ = 1249'K 2/3 'W 4/3 So /3 A1 1/6 JrrTV 1 . (107) 

The turbulent velocity of the coronal plasma may be found from 
relation 



(108) 



Only the perpendicular fluctuating velocity is taken into ac- 
count since we are only concerned with shear- Alfven waves in a 
strongly anisotropic turbulence. We substitute P x by its expres- 
sion d73l to finally obtain the following prediction: 



u± =50W 1/2 <Z/ 1/3 £ 2 / 3 AT 1/6 £- 1/2 kms- 



(109) 



These results compare favorably with observations, in partic- 
ular, in the quiet corona with a heating flux large enough to ex- 
plain the observations. Note that the heating prediction is mostly 
sensitive to the magnetic field intensity So which has the larger 
power law index by 5/3. Thus, S = 3-4 (B = 3-4 x 10~ 2 G) 
leads to a factor about ten times larger for the heating flux with a 
value close to the measurements 10 4 Jm _2 s _1 for active regions. 

5. Coronal hole heating 

5. 1 . Geometry of open magnetic field lines 

In this Section, the model prediction for the heating rate is ex- 
tended to coronal holes where the plasma is guided along a 
large-scale magnetic field which expands into the interplanetary 
medium. It is along such structures (mainly at the poles) that 
the fast solar wind is released whereas the slow wind is freed at 
lower latitudes around the equatorial plane. In these configura- 
tions of open magnetic lines, the reflection of outward Alfven 
waves, due to some inhomogeneities, produces inward waves 
which eventually sustain nonlinear interactions. In this case, 
the coronal heating is clearly dependent on the reflection rate 
of Alfven waves jVelli et al. l"993llDmitruk & Matthaeus 2003l >. 
The precise origin of the partial reflection of Alfven waves is 
still under debate. Therefore, any theoretical prediction is use- 
ful for the comparison between observations and models, and 
eventually for the understanding of the solar corona dynamics. 
In this Section, we seek theoretical predictions like the reflec- 
tion rate needed to sustain an efficient coronal heating localized 
in coronal holes. We see that a small reflection rate of 15 - 20% 
is enough to recover the coronal heating observations. This re- 
sult may constrain the efficiency of the mechanisms invoked to 
produce reflected Alfven waves. 

Figure|5]shows a schematic view of open magnetic field lines 
which only differs from the magnetic loop configuration (Figure 
|3j is the upper boundary condition. This second photospheric 
surface is replaced by a permeable boundary from which out- 
ward Alfven waves may be partially reflected. Note that the first 
model of magnetic arcades corresponds to the case where out- 
ward Alfven waves are totally reflected and where outward and 
inward waves are balanced. 
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Fig. 5. Coronal hole geometry: the coronal plasma is confined in 
a volume delimited by —h < x < +h and -oo < y < +00, and may 
escape from the low coronal medium at the top of the structure 
(for z > £)■ The photosphere appears as a (single) lower bound- 
ary plane at altitude z = The interplanetary inhomogene- 
ity allows the reflection of a fraction of outward Alfven waves 
which produces inward waves. 

5.2. Reflection rate versus cross-helicity 

To quantify the reflection rate of Alfven waves (also called the 
Alfvenicity), a characteristic quantity is the normalized cross- 
helicity cr c defined as 

z +2 -z- 2 
z +z +z~ z 

where z + is the Alfven inward wave and z~ is the outward wave 
(with z s = |z s |, s = +). Such a quantity provides a measure of 
the relative importance between outgoing and ingoing Alfven 
waves. In particular, cr c = +1 involves the presence of only one 
type of polarity, and consequently, the absence of nonlinear in- 
teractions (thus no heating at dissipative scales), whereas a bal- 
ance between z + and z~ waves implies cr c — leading to more 
important nonlinear interactions (between waves of different po- 
larity) as was assumed in the magnetic arcade (loop) configura- 
tion analyzed in previous Sections. The reduced cross-helicity is 
indeed related to the wave reflection rate, %, such as 



The case "R — 1 (all outward waves are reflected) is similar to the 
loop configuration in which there is the same number of inward 
and outward waves, whereas in the case of 7? — no waves are 
reflected. 

5.3. Unbalanced turbulence for coronal holes 

The integro-differential kinetic equations d53l l for the Alfven 
waves are numerically integrated with a logarithmic subdivision 
of the k axis : kj = 8k2 l ' F , where i is an integer (i = 1,N), 
6k = 0.125 is the smallest wave number reached in the compu- 
tation, and F = 8 is the refinement of the grid. (This method 
was previously performed and optimized in Ga ltier et al. 20001 ) 
This numerical technique allows to reach larger Reynolds num- 
bers since, for a wave number resolution N = 157, the maxi- 
mum wave number reached, k lmx , is around 10 s (the magnetic 



Prandtl number is equal to unity). Note that dissipative terms 
are introduced in the inviscid kinetic equations to avoid numer- 
ical instabilities, and the viscosity is fixed to 10~ 5 . The initial 
Elsasser fields are injected in the wave number range [<5&,44] 
with an energy spectrum proportional to k\ exp(— peak- 
ing at k ±0 ~ 1. The corresponding U rms and integral scale Lq are, 
respectively, about 4.75 and 1.5, which gives an initial Reynolds 
number (U rm sLo/v) of about 10 6 . The flow is then left to freely 
evolve. 

In Figure |6l we show instantaneous energy spectra of the 
shear- Alfven waves, and E~, obtained for, respectively, z+ 
and z~]_ fluctuations with a cross-helicity of 0.8 (7? = 25%) (In 
the following we respectively consider z~]_ and z\_ as the inward 
(reflected) and outward waves). The chosen time is the one at 




10 10 1 10 2 10 3 10 4 10 5 



Fig. 6. Energy spectra of shear- Alfven waves E\_ (solid line) and 
E~ ± (dashed line) with % = 25% (cr c = 0.8). The straight line 
follows a k] 2 law which is the theoretical prediction for the bal- 
anced turbulence case (7? = 100% and cr c — 0). Inset: compen- 
sated product of energy spectra as EJ_ k\. 



which the energy spectra are the most developed, ie. the most 
extended towards large wave numbers k±. Different behaviors 
are clearly found for the E]_ and E]_ spectra with, at large scales, 
a domination of inward Alfven waves linked to the choice of the 
initial conditions. At smaller scales, we see the appearance of 
an extended inertial range where the inward waves (£") slightly 
dominate. Note that the energy spectra product, E]_El_, follows 
a kj 4 scaling law as predicted theoretically (G altier et al. 2000"i l. 
and shown in the inset of Figure|6] Finally the two energy spectra 
overlap at dissipative scales. From such a numerical simulation, 
we are able to compute the turbulent viscosities for unbalanced 
turbulence by integration of the spectra over wave numbers, from 
the beginning of the inertial range up to the dissipative scales. 
This computation is done for several simulations corresponding 
to different reflection rates. 
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5.4. Parametric study and predictions 

A parametric study of the heating flux according to the reflection 
rate of outward waves is performed. The turbulent viscosities are 
calculated for different values of the reflection rate by integrating 
the shear- Alfven wave spectra, from the injection scale (normal- 
ized to unity in the numerical simulation) up to the dissipative 
scales, and by using relation d62b to estimate v s ± . In such a calcu- 
lation, the chosen time is the one at which the energy spectra are 
most developed. Since the result depends on the initial amount 
of energy taken in the simulation, the flux is normalized to the 
flux obtained in the loop case (see Section |4| for which a bal- 
anced turbulence is assumed. With this method, we are able to 
approximately predict the variation of the relative heating flux as 
a function of the reflection rate. 

From the turbulent flux expression d75l l, we obtain the rela- 
tions 

\T:\<R)\ ( v s \W) \ 2 _ [jEl s CR)d K± ] 2 

\n \ ^ I = [fE±dK± f ' (112) 

which depend on the reflection rate of outward waves %. T- and 
T z are, respectively, the heating fluxes in the unbalanced and 
balanced turbulent cases. The presence of the directional polar- 
ity s — ± keeps track of the advection between waves of dif- 
ferent polarity. The calculation made here concerns the heating 
flux at the lower boundary (z = -€) since the upper boundary 
condition previously used is no longer the same. We have also 
assumed, in particular, that the Kolmogorov constant does not 
change drastically for different reflection rates. Figure [7] shows 
the flux ratio as a function of the reflection rate %. As expected, 
the curve decreases with the reflection rate, meaning that the tur- 
bulent heating is clearly less efficient for smaller values of %. 
Finally, the flux ratio goes to zero when only one type of wave 
remains (since in this case the nonlinear interactions disappear). 
Note that the decreasing function is close to a linear curve with 
a flux ratio of about 40% for % = 50%. 

The heating flux prediction in coronal holes is known and 
estimated as around 100Jm _2 s _1 ( Withbroe et al. 19771 >. Since 
the background magnetic field is still about 10~ 2 T, we may find 
the reflection rate in coronal holes by using both the relation 
( 1107b and the curve in Figure [7] The value % — 25% leads to a 
flux ratio of about 15%— 10%, a small enough value to recover 
the heating flux for coronal holes (without taking into account 
the fast solar wind). We see that a relatively small reflection rate 
is therefore enough to produce an efficient coronal heating due 
to anisotropic MHD turbulence. 

6. Discussion and conclusion 

In this paper, we have developed an analytic model for strongly 
anisotropic structures in order to recover well known coronal 
heating rates for the quiet sun, active regions and coronal holes. 
The coronal structures are assumed to be in a turbulent state 
maintained by the slow motion of the magnetic footpoints an- 
chored on the photospheric surface. The main difficulty is that 
existing spacecraft are unable to resolve all the inertial (small) 
scales, and a fortiori the dissipative scales. So, firstly, the in- 
jected energy at large scale is resolved from incompressible 
MHD equations, and for specific boundary conditions satisfying 
the divergence free condition. Secondly, the (unresolved) small- 
scale dynamics is modeled by turbulent viscosities derived from 
an asymptotic (exact) closure model of wave turbulence for the 
case of the loop configuration (also called balanced turbulence), 
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Fig. 7. Evolution of the heating flux according to the reflection 
rate of outward Alfven waves. The heating flux is normalized to 
the flux obtained in the balanced turbulence case (Section|4|, i.e. 
for?? = 100%. 



i.e. when there are as many inward as outward waves which non- 
linearly interact. In the open magnetic line configuration the non- 
linear interactions are sustained by the conversion of outward 
waves into inward waves by reflection (unbalanced turbulence). 
We have numerically integrated the kinetic equation of Alfven 
waves for different values of the reflection rate, and plotted the 
evolution of the coronal heating flux according to this rate, as 
a measure of the relative magnitude of ingoing and outgoing 
waves. Finally these results are compared to the balanced tur- 
bulence flux (loop configuration). 

For standard loop geometry parameters, and for a magnetic 
field intensity of Bo = 10~ 2 T, we find a heating flux prediction 
of about 1.2 x 10 3 J m~ 2 s _1 which is close to the value measured 
in the quiet sun (10 3 Jm _2 s _1 ). Moreover, the prediction for the 
turbulent velocity (50 km s _1 ) compares favorably with the mea- 
surements of nonthermal velocities (30 km s ) in the quiet solar 
corona, as well as with the maximal values of some line profiles 
(55 km s" 1 ) found by the SUMER instrument (IChae et al. 19981 >. 
For active regions, our estimation tends toward a value of about 
10 4 Jm _2 s _I for a magnetic field intensity stronger (3-4- times 
larger) than for the quiet sun, which is an expected value for 
strong solar activity. For the heating rate prediction in open mag- 
netic lines, we take as a reference the loop configuration (for 
which there is a balanced turbulence) and the same parameters 
as for the quiet sun. We can recover the well known estimates 
for coronal holes of about 100 Jm^V (without taking into ac- 
count the fast solar wind) for a reflection rate of outward waves 
of about 15-20%. 

The previous predictions do not take into account the coro- 
nal heating by pseudo-Alfven waves (nor the heating due to the 
pure 2D state - with wave vectors such as (k.±,k\\ — 0) - which 
is not described by Alfven wave turbulence), since their par- 
allel fluctuations are excluded by the boundary conditions (in- 
volving only perpendicular fluctuations). Thus, their insertion 
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should heat slightly more the solar corona and raise our predic- 
tions slightly. Moreover, in the computation of heating flux for 
the unbalanced turbulence case, only the heating due to outward 
Alfven waves is considered. The heating due to inward waves 
should modify slightly the heating prediction for coronal holes, 
and therefore allows a slight decrease of the reflection rate with 
a heating rate maintained around at 100 J m _2 s _1 . 
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